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1 .  Introduction 


Fiber  reinforced  composites  have  been  used  for  about  two  decades  for  improving  the  burst  (due 
to  internal  pressure)  and  collapse  pressure  (due  to  external  pressure)  of  spherical  pressure  vessels 
[1-3].  Although  several  reinforcing  patterns  had  been  tried  in  the  past,  for  example  radial 
filament  reinforcement  [3],  the  accepted  technology  for  reinforcing  spherical  vessels  has  emerged 
to  be  the  quasi-isotropic  lay-ups  in  the  0-(p  plane  (0:  colatitude,  <j):  azimuth)  (Figure  1). 


x 

Figure  1 .  Coordinate  System 

Gerstle  [4|  has  reported  the  stress  and  strength  analysis  of  single-layered  fiber  reinforced 
spherical  vessels  with  elastic-ideally  plastic  bladder.  In  his  work  the  layer  of  a  vessel  is 
considered  to  be  composed  of  quasi-isotropic  lay-ups  of  one  type  of  material.  Gerstle  has  found 
that  the  burst  pressure  of  a  single-layered  vessel  is  limited  by  the  radial  compressive  strength  of 
the  vessel.  The  radial  compressive  strength  of  these  vessels,  incidentally,  is  related  to  the 
stiffness  and  transverse  compressive  strength  of  the  laminae  used  in  the  lay-ups.  Thus  one 
logical  way  of  improving  the  failure  pressure  is  to  make  hybrid  vessels  where  an  optimal  use  of 
the  stiffnesses  and  strengths  of  the  layers  is  expected  to  improve  the  failure  pressure.  Gerstle  [5] 
has  then  presented  a  similar  work  for  two-layered  hybrid  spherical  vessels  and  has  shown  a  way 
to  use  the  stiffnesses  and  strengths  of  the  two  hybrid  layers  in  improving  the  burst  pressure  over 
that  of  a  single-layered  vessel.  The  work  presented  here  is  in  a  way  an  extension  of  Gerstle's 
work  and  is  applicable  to  any  number  of  layers. 


2 .  Equations  for  Stresses  and  Strains  for  Multilayer  Spheres 


The  effective  properties  of  the  layers  of  a  spherical  vessel  with  quasi-isotropic  lay-ups  are 
transversely  isotropic  with  isotropy  in  the  the  &-(j)  plane  (Figure  1).  The  vessel  is  considered 
spherically  symmetric.  The  stress- strain  relations  of  the  i-th  layer  of  the  vessel  reduce  to  [6] 


<*°  =  C<V>  + 

(1) 

(2) 

=  +  of 

(3) 

(4) 

P(«)  _  durl)  .  p(0  _  Url)  .  p(0  _  JO 
*  dr  ’  «  “  r  ’t(P 

(5) 

Here  Cifl\  S^K  ( i,j  =  r,  6,  <p)  are  the  elements  of  the  three-dimensional  effective  stiffness  and 
compliance  matrices  of  the  i-th  layer  respectively.  A  method  of  determining  the  three- 
dimensional  effective  stiffnesses  and  compliances  of  a  laminate  is  given  in  reference  [7]  and  is 
summarized  in  Appendix  A.  If  the  sphere  is  subjected  to  both  internal  and  external  pressure, 
then  the  boundary  conditions  are 

when  r  =  a  or  V)  =  -p  (the  negative  sign  is  because  of  the  compressive  pressure)  (6) 

when  r  =  b  op )  =  -q  (the  negative  sign  is  because  of  the  compressive  pressure)  (7) 

At  the  interface  between  two  adjacent  layers  the  following  conditions  must  be  satisfied 

at  the  interface  r  =  a,  <r/d  =  or  <i+l)  (8) 

up  =  ur<i+I>  (9) 

If  we  express  the  stresses,  strains,  and  displacement  in  terms  of  the  interlayer  normal  tractions, 
qh  (Figure  2) 
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Figure  2.  Cross  section  of  the  sphere  and  interlayer  tractions,  qj,  used  in  the  analysis, 
then  the  expressions  for  stresses,  strains,  and  displacement  of  the  sphere  are: 


a(l)  =  -1 
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(15) 


e(0  =  e« 
ce 


Qi-ic*+*-qi 

l-c-2”' 

(n,-i)C®+2C® 

where 


<3.  i  1 

Ci  =  iT’  n‘  2 


r(i)  r(0 


Incidentally,  when  r  =  a  equation  (10)  becomes 

(r  =  aQ  =  a)  =  -qQ  =  -p  (inside  pressure)  (17 

and  when  r  =  b  equation  (10)  again  reduces  to 

ojn)  (r  =  an  =  b)  =  -qn  =  -q  (outside  pressure)  ( 1 8 

By  expressing  the  stresses  in  terms  of  the  interlayer  normal  tractions,  <fr,  equation  (8)  is 
automatically  satisfied.  Then  substituting  equation  (16)  into  equation  (9)  we  get  a  set  of  linear 
equations  to  determine  the  interlayer  normal  tractions,  q{. 

Ql+l  r‘+1  +  gi(A+i-A-K+i-K)g«  +  qi- 1  ociOi-i  =  o  (19 

i=  1,2,  .  .  .  ,  (n-l) 

here,  q0=p  (internal  pressure),  qn  =  q  (external  pressure)  are  known,  and 
2 


r'  “  /li 
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The  interlayer  normal  tractions,  {qi,  i  =  1 , .  .  . ,  n-1),  are  known  by  solving  equation  (19).  Then 
substituting  the  values  of  qL  in  equations  (10-16),  the  stresses,  strains,  and  displacement  for  all 
layers  are  determined. 

3 .  Stresses  and  Strains  for  Single  Layer  Spheres 

If  there  is  only  one  material  for  the  spherical  shell,  super  and  subscripts  in  the  equations  above 
assume  only  one  value;  i.e.,  i  =  1;  then  in  the  inside  radius,  r  -  a0-  a\  and  on  the  outside  radius, 
r  =  a\  -  b_  We  can  define  the  ratios  of  the  radii  as: 


All  equations  in  the  previous  section  are  simplified  for  the  case  of  one  layer  of  material  which  can 
also  represent  a  homogenized  laminated  shell. 

4 .  Failure  Criterion  for  Three-Dimensional  Laminates 

Failure  criteria  for  a  three-dimensional  body  can  be  easily  written  down  in  terms  of  the  most 
general  quadratic  form.  There  will  be  a  maximum  of  15  strength  parameters. 

A  transversely  isotropic  materia*  is  a  reasonable  representation  of  a  unidirectional  composite  ply  if 
the  diameter  of  the  fiber  is  small  relative  to  the  thickness  of  the  ply  and  there  are  many  fibers 
randomly  distributed  in  the  ply.  With  this  symmetry,  the  number  of  strength  parameters  can  be 
reduced  to  eight.  If  we  can  further  invoke  the  identity  between  tension-compression  stresses  and 
pure  shear  oriented  45-degrees  from  the  former  stresses,  this  additional  relation  further  reduces  the 
strength  parameters  in  a  3-D  transversely  isotropic  materials  to  seven. 
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•  FjjCTiaj+FiCT^  =  1,  for  3-D  contracted  notation:  crq  =  cyz,  ar  =  azx,  os  =  axy 


Fxxax+Fyya5+Fzzo?  +  Fqqaq+Frrcj?  +  Fsso| 


No.  parameters- 


+  2Fxyaxay  +  2Fxzaxaz+2Fyzaycjz+Fxax+Fyay+Fzaz+Fqaq+FrCTr+FsCTs  =  1  — 

From  symmetry:  Fyy  =  Fzz,  Frr  =  Fss,  Fxy  =  FXZ/  Fy  =  Fz,  Fq  =  Fr  =  Fs  =  0 - 

Fxx^x  +Fyy(ay  +  az  )+Fqqaq +  Fss(Or +cj|) 

+  2Fxy<jx(ay+az)+2Fyzayaz+Fxax+Fy(ay+az)  =  1 


From  transverse  isotropy:  oy  =  -az  =  Q,  other  <jj  =  0 
Fyy(2Q2)-2FyzQ2  =  1 ,  or  Q2  =  -  or  Q  = 
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vw 


2(Fyy-Fyz) 


72C1-F*) 


where  Fyz=  Fyz s/ Fy yFzz  -  FyzFyy,  Fyz  =  Fyz/Fyy 


Figure  3.  3-D  quadratic  failure  criterion  for  a  transversely  isotropic  material 

Thus  the  difference  between  2-  and  3-D  quadratic  failure  criteria  for  a  unidirectional  ply  is  by  only 
one  additional  strength  parameter.  The  reduction  in  the  number  of  strength  parameters  from  15  to 
seven  is  illustrated  in  Figure  3. 

5 .  Strength  Ar  lysis  of  Spheres 

Stress  analysis  of  a  spherical  shell  in  3-D  is  similar  to  that  for  2-D  provided  appropriate  elastic 
moduli  and  stress-strain  relations  are  used.  3-D  analysis  will  of  course  require  3-D  moduli  and  3- 
D  stress-strain  relations.  This  is  illustrated  in  Figure  4  where  over  bars  represent  the  average  stress 
or  strain  of  the  3-D  thick  laminate.  The  relation  between  laminate  stress  and  ply  stress  is 
determined  with  the  knowledge  of  the  laminate  compliance  and  ply  stiffness,  both  in  3-D. 


{e}  =  [S]{a},  where  [S]  is  the  3-D  laminate  compliance 

{a}  =  [C]{e),  where  [C]  is  the  3-D  ply  stiffness 

lFijCTi<Jj]R2+[Fjaj]R-  1  =  0, 

where  F*j,  F,  are  3-D  strength  parameters 
Assuming  unidirectional  plies  are  transversely  isotropic: 

Fxx^x  "*’Fyy(0y  +  2FyzGy0z'*Cr2  +  30q)  +  2F xy^X^^y+^z) 

+  Fss(a^  +  a?)+Fxox+Fy(CTy+a2)  =  1 

Figure  4.  Relations  between  ply  stress  and  laminate  stress  in  a  3-D  thick  laminate.  The 
quadratic  failure  in  3-D  can  be  applied  on  a  ply-by-ply  basis. 

In  particular,  the  uniaxial  tensile  or  compressive  stress  applied  along  the  z-axis  can  be  easily 
calculated,  and  the  resulting  failure  strength  can  also  be  calculated  using  the  quadratic  criterion. 
The  interaction  term  in  the  isotropic  plane  F*yz  can  be  show;,  to  have  a  strong  influence  on  the 
uniaxial  compressive  strength.  The  back-calculated  interaction  term  for  a  typical  graphite-epoxy 
composite  is  approximately  -0.84.  This  back-calculation  is  shown  in  Figure  5  for  both  thick  cross- 
ply  and  (±45]  laminates. 


-1.0  -0.9  -0.8  -0.7  -0.6  -0°5  -1.0  -0  9  -0.8  -0.7  -0.6  -0°5 


Normalized  interaction  term,  F*z 

Figure  5.  The  best-fit  interaction  term  in  the  isotropic  plane  based  on  the  measured 
uniaxial  compressive  strength. 


It  is  seen  that  the  quadratic  failure  criterion  required  a  highly  coupled  relation.  The  maximum 
stress  or  maximum  strain  criterion  assumed  that  interaction  between  failure  modes  does  not  exist. 
Based  on  the  test  results  shown  in  Figure  5,  the  predicted  maximum  compressive  strength  would 
be  the  same  as  the  transverse  compressive  strength  in  a  unidirectional  ply.  The  value  for  most 
graphite-epoxy  composites  would  be  in  the  range  of  200  MPa.  Thus  an  error  of  300  percent  is  the 
result  of  ignoring  failure  mode  interactions.  We  therefore  believe  that  the  quadratic  criterion  is 
flexible  in  taking  into  account  such  mode  interactions;  while  the  maximum  stress  or  maximum 
strain  criterion  is  rigidly  defined  and  does  not  recognize  the  interactions. 

6 .  Burst  and  Collapse  Pressure  of  Spheres  Made  of  One  Material 

The  burst  and  collapse  pressures  of  a  thick  spherical  shell  under  internal  and  external  pressures 
based  on  strength  can  be  derived  using  the  stress  analysis  for  thick  shells  coupled  with  appropriate 
failure  criteria.  In  the  present  study,  it  is  assumed  that  the  shell  failure  is  controlled  by  strength, not 
by  buckling.  It  is  further  assumed  that  the  shell  consists  of  homogeneous,  quasi-isotropic  material 
with  3-D  elastic  constants  listed  in  Appendix  A. 

The  predicted  burst  pressures  for  typical  composite  materials  are  shown  in  Figure  6.  From  the 
equations  of  stress  and  strain,  it  is  observed  that  they  are  dependent  on  the  ratio  of  the  outside  to 
inside  radii,  not  on  their  absolute  values.  The  displacements  on  the  other  hand  are  functions  of  the 
absolute  radii. 
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Burst  presssure,  MPa 


Figure  6.  The  burst  pressures  of  two  graphite-epoxy  composite  spheres  as  functions  of 
the  ratios  of  radii. 

Open  symbols  for  each  material  are  shown  for  failure  on  the  outside  diameter  and  solid  symbols  on 
the  inside.  For  thin  shells  failure  occurs  on  the  outside;  for  thick  shells  failure  occurs  on  the 
inside.  Both  materials  reach  some  asymptotic  strength  level,  but  their  functional  dependences  on 
the  ratio  of  radii  are  different.  Below  a  ratio  of  1.35,  IM6  has  higher  burst  pressure;  above  this 
ratio,  T300  becomes  higher.  It  is  therefore  not  possible  to  guess  what  material  would  be  stronger 
unless  a  calculation  is  made  First. 

7.  Burst  Pressure  of  Hybrid  Spheres  Made  of  Two  Materials 

« 

One  of  the  objectives  of  this  work  is  to  design  a  composite  sphere  to  operate  at  200  MPa.  For  the 
initial  design  a  safety  factor  of  1.5  was  considered.  Thus  the  design  of  a  sphere  that  could 
withstand  a  maximum  pressure  of  300  MPa  was  considered  to  be  a  safe  design.  It  was  initially 
inferred  from  Figure  6  that  a  sphere  made  with  T300/5208  composite  of  b/a=1.5  could  take  the 
pressure  close  to  300  MPa.  As  a  part  of  a  demonstration  of  the  design,  a  T300/epoxy  sphere  of 
inner  diameter  of  240  mm  and  b/a=l  .5  was  constructed.  The  sphere,  however,  failed  prematurely. 
A  close  examination  of  the  video  photographs  of  the  test  indicated  that  the  failure  initiated  at  the 
interface  between  the  liner  nozzle  and  the  composites.  For  this  sphere  of  b/a=1.5,  the  thickness  of 
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the  composite  was  approximately  60  mm.  Relating  to  the  complex  winding  pattern  involved  in  the 
vicinity  of  the  liner  nozzle  of  the  sphere,  the  manufacturer  (Courtaulds  Advanced  Materials, 
France)  pointed  out  the  difficulty  in  maintaining  a  good  structural  integrity  of  the  interface  between 
the  liner  nozzle  and  the  composites  for  such  a  thickness  of  the  sphere.  The  above  premature  failure 
may  thus  be  attributed  to  the  poor  interface  between  the  nozzle  and  composites.  Incidentally,  two 
IM6/epoxy  spheres  of  b/a=1.25  and  1.33  of  inner  diameter  of  240  mm  were  successfully  tested 
earlier  by  CEA  (Commissariat  a  l'Energie  Atomique,  France)  where  failure  pressures  were 
respectively  within  10  and  4%  of  the  predicted  pressure.,  as  shown  in  Figure  7.  The  failure 
pressures  of  these  spheres,  incidentally,  were  much  less  than  our  present  design  objective.  One 
possible  way  of  Increasing  the  failure  pressure  keeping  the  same  b/a  ratio  is  to  consider  the  design 
of  hybrid  spheres.  Thus  in  view  of  the  previous  success  in  manufacturing  sphere  of  b/a=1.25,  a 
design  restriction  of  limiting  b/a=1.25  was  imposed  for  the  design  of  hybrid  spheres. 


1.0  1.1  1.2  1.3  1.4  1.5 

OUTER  RADIUS 
INNER  RADIUS 

Figure  7.  Comparison  of  the  predicted  burst  pressure  with  the  tested  data.  Material:  IM6/Epoxy 

The  stress  and  strength  analysis  discussed  above  is  applied  for  the  design  of  two-layered  hybrid 
spheres  made  of  T300/5208  and  IM6/epoxy  composites.  Both  the  layers  are  of  quasi-isotropic 
laminates.  The  predictions  of  the  burst  pressure  as  functions  of  percent  of  layer  thickness,  m,  are 
shown  in  Figure  8.  The  notation  [(T300)m/(IM6)nl  indicates  that  the  inner  layer  is  made  of 
T300/5208  composite,  and  thickness  of  the  inner  layer  is  m  percent  of  the  total  thickness  of  the 
sphere.  Obviously,  the  thickness  of  the  outer  layer  is  ( 100-m )  percent  of  the  total  thickness.  Thus 
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the  sum  of  m  and  n  becomes  100.  When  either  m  or  n  becomes  zero  or  100,  the  sphere  is  no 
longer  a  hybrid  sphere  but  rather  a  sphere  of  one  material.  It  is  obvious  from  the  figure  that  a 
hybrid  sphere  with  T300  inside  and  IM6  outside  yields  a  higher  burst  pressure  than  a  sphere  with 
reversed  order  of  the  layers.  It  is  also  noticed  from  the  figure  that,  for  the  hybrid  sphere  with  T300 
inside,  if  the  thickness  of  either  of  the  layers  is  varied  between  20  and  80%  of  the  total  thickness, 
the  variation  of  the  burst  pressure  is  limited  to  within  only  2.5%  of  240  MPa.  Thus  within  the 
range  of  20  to  80%  of  the  layer  thickness,  any  small  manufacturing  error  in  the  layer  thickness  will 
hardly  alter  the  expected  burst  pressure. 

b/a=1.25,  Quasi-isotropic 


LAYER  THICKNESS  PERCENT,  m 

Figure  8.  The  burst  pressure  of  two  hybrid  spheres  as  functions  of  the  layer  thickness. 

8.  Conclusions 

We  have  presented  a  procedure  to  determine  the  stress,  strain,  and  ultimate  internal  and  external 
pressures  of  a  homogenized,  isotropic  spherical  shell  of  any  thickness.  It  can  be  shown  that  the 
membrane  solution  would  correspond  to  the  tangent  of  the  burst  pressure  curves  in  Figure  6 
above. 

Comparable  curves  for  other  materials  will  require  the  3-D  elastic  moduli  and  seven  strength 
parameters  needed  for  the  quadratic  failure  criterion.  The  collapse  pressure  of  spherical  shells 
under  external  pressure  can  be  determined  precisely  the  same  way.  Only  boundary  conditions  need 
to  be  changed. 
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The  required  calculation  can  be  easily  programmed  in  a  small  computer.  In  fact,  a  Fortran  program 
is  written  to  perform  the  numerical  analysis  for  multilayered  (or  hybrid)  spheres.  The  instructions 
for  executing  this  code  are  given  in  appendix  B. 
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Appendix  A 


Method  of  Predicting  Three-Dimensional  Effective  Moduli  of  Laminates: 

The  sphere  is  considered  to  be  made  of  layers  with  quasi-isotropic  lay-ups.  The  stress  analysis  is 
performed  by  considering  the  three-dimensional  effective  moduli  of  each  layer.  The  method  of 
predicting  an  equivalent  homogeneous  effective  moduli  of  a  laminate  is  discussed  in  detail  in 
reference  [7].  Briefly,  the  method  is  based  on  matching  the  boundary  displacements  of  the 
homogeneous  system  with  those  of  the  laminated  system  of  three  interrelated  boundary  value 
problems  (BVP).  The  BVP  are:  (a)  ring  subjected  to  equal  hydrostatic  pressure  on  both 
surfaces,  (b)  cantilever  beam  under  end  transverse  load,  and  (c)  cantilever  beam  subjected  to 
uniformly  distributed  load  on  top  and  bottom  surfaces.  These  interrelated  BVP  are  selected  to 
predict  the  effective  moduli  for  in-plane  and  bending  mode  of  deformation.  The  analyses  for 
predicting  the  effective  moduli  are  given  below: 

Problem  1 .  Ring  with  Hydrostatic  Pressure 

This  problem  is  used  to  predict  the  interlaminar  effective  moduli  of  symmetric  and  orthotropic 
laminates  for  in-plane  extension.  For  a  large  value  of  the  radius  to  thickness  ratio,  r/t  (r:  radius, 
t:  thickness),  the  deformation  of  the  ring  approximately  represents  the  in-plane  deformation.  The 
laminated  ring  and  its  equivalent  homogeneous  ring  are  subjected  to  equal  hydrostatic  pressure 
on  both  inner  and  outer  surface,  p,  as  shown  in  Figure  9.  The  inner  and  outer  displacements  of 
the  homogeneous  ring  are  matched  with  those  of  the  laminated  ring  to  obtain  effective  moduli  of 
the  laminate  of  the  ring.  The  interlaminar  effective  elastic  constants  obtained  from  this  problem 

_0  -o  -o 

are  ^3,  v3 1,  and  v32,  where  Vjj  is  the  Poisson's  ratio  for  Ei  under  oj.  The  superscript  'o’  is  used 
to  indicate  that  the  elastic  constants  are  associated  with  the  in-plane  deformation. 


14 


b 


— *1 

(a)  LAMINATED  RING 


(b) 


I - b — H 

EQUIVALENT 
HOMOGENEOUS  RING 


Figure  9.  Configuration  of  the  laminated  and  its  equivalent  homogeneous  ring  subjected  to  equal 
hydrostatic  pressure  on  the  inner  and  outer  surface. 


The  method  of  obtaining  the  radial  displacement  of  the  laminated  composite  ring,  u'(r), 
(Figure  9a)  is  given  by  Lekhnitskii  (1968).  The  subscript  'i'  indicates  the  j-th  layer  of  the 
laminated  ring.  We  do  not  find  the  necessity  of  repeating  the  procedure  for  obtaining  the 
expression  of  u'(r)  in  this  report.  The  expression  for  the  radial  displacement  of  the  equivalent 
homogeneous  ring  (whose  effective  moduli  are  to  be  determine),  Figure  9b)  is  similar  to  that  of 
any  layer  of  the  laminated  ring;  that  is: 


(Al) 


where  m  =  £■,  k  =  a:  inner  radius,  b:  outer  radius,  and  S;;  (i,j=2,3)  are  the  elements 

b  V  S22 

of  the  effective  compliance  matrix. 

By  equating  the  inner  and  outer  displacements  of  the  homogeneous  ring  with  those  of  the 
laminated  ring,  we  get  two  nonlinear  equations  in  k,  S22,  and  S23  as  follows: 


-~r  fS22k(l+m2k-2mkl)-S23(l-m2k)]-u1(a)  =  0  . 

l-m2kL  J  ,  (forr=a,  t=l)  (A2) 


7— S22k(2mk+1-m2k-l)-S23(l-m2k)l-un(b)  =  0 

l-m2kL  J  ,  (for  r=b,  i=n) 


(A3) 


The  laminated  ring  is  assumed  to  have  n  number  of  layers  in  equation  (A3).  Here  S22  is 
considered  to  be  known  from  the  laminated  plate  theory  which  seems  to  be  a  reasonable 
assumption.  Then  S23,  k,  and  hence  S33  are  solved  from  equations  (A2)  and  (A3)  by  the 
Newton-Raphson's  technique.  The  laminate  effective  properties  in  terms  of  S22,  S23>  and  S33 
arc: 

p°_  J_  ~°  S23 

3  S„-  v«-  S22  (A4) 

— o 

The  constant  v3i  can  similarly  be  obtained  by  rotating  the  laminate  by  90°  about  the  r-axis  of  the 
ring  (a  90  degree  rotation  brings  the  1-axis  of  the  laminate  with  the  0-direction  of  the  ring).  The 
superscript  'o'  is,  again,  to  denote  that  the  effective  properties  apply  for  in-plane  deformation. 

Problem  2:  Beam  with  Transverse  End  Load 

The  laminated  beam  is  subjected  to  an  end  transverse  load,  P,  which  is  modeled  as  an  end 
transverse  shear  load.  The  displacement  field  of  the  laminated  beam  is  obtained  by  considering 
the  beam  in  the  state  of  plane  stress  either  in  the  1-3  or  in  the  2-3  plane.  The  plane  of  lamination, 

rf  rf 

however,  is  the  1-2  plane  in  both  cases.  From  this  problem  the  effective  properties,  u  13,  u23, 

-f  _f 

v3\,  and  v32  are  estimated.  The  superscript  'f  is  used  to  indicate  that  the  properties  are 
associated  with  bending  as  the  dominant  mode  of  deformation. 

LOAD, P  LOAD, P 


(c)  ENLARGED  VIEW  OF  BOUNDARY  CONDITIONS  AT  POINT  A 

Figure  10.  Configuration  of  the  laminated  beam  and  its  equivalent  homogeneous  beam  under 
end  transverse  load. 

The  expressions  for  the  displacement  field  can  easily  be  obtained  by  starting  with  an  Airy 
stress  function  appropriate  to  this  problem.  Since  the  end  transverse  load,  P,  is  modeled  as  the 
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transverse  shear  load,  (see  Figure  10),  it  is  logical  to  assume  that  az  is  zero  everywhere  in  the 
beam.  The  Airy  stress  function,  Fj(x,z),  for  a  constituent  orthotropic  lamina  of  the  above 
problem  for  stresses  in  the  (x-z  or  1-3)  plane,  has  the  form 


Fj(x,z)  =  Ajxz  +  BjZ2  +  Cjxz2  +Djz3  +  Ejxz3  (A5) 

where  the  constants,  Aj,  Bj, . . .,  Ej  are  associated  with  the  properties  of  the  j-th  lamina.  After 
using  the  linear  stress-strain  constitutive  relations  and  the  expressions  for  stresses  in  terms  of  the 
above  stress  function,  Fj(x,z),  the  displacement  fields  of  the  laminated  beam,  uj(x’z)  and 
wj(x,z\  along  the  x  and  z  directions  respectively,  are: 


»?<*.*>  =  f +(S®+S“>)z2]+Cj[-S®x2z+(S«+S«)z3] 


+DJS(j’1)xz+EJS^jx+QJz  +  RJ 


(A6) 


wj (x,z)  =  AjSjjX - B jS^xz + Cj -t-S^xz2  J 
+iDj[snz2-sn  x2]+EjS^3 z-Qjx+Mj 


CAT) 


The  constants  Aj,  Bj, .  .  .  ,  Ej,  Mj,  Qj,  Rj,  for  the  j-th  lamina,  are  determined  from  the  contact 
conditions  at  the  interface  of  each  lamina,  boundary  traction,  P,  and  the  boundary  conditions  at 
the  support  end,  x=L.  Here  Sij(k)  (i,j  =  1,3,5)  are  the  elements  of  the  compliance  matrix  of  the 
k-th  lamina. 

Following  a  similar  procedure  for  equations  (A6)  and  (A7),  the  expressions  for  the 
displacement  fields,  u^(x,z)  and  w^(x,z),  for  the  equivalent  homogeneous  beam  (Figure  10b)  for 
plane  stress  in  (1-3)  or  (x-z)  plane,  keeping  a  horizontal  element  at  point  A  (Figure  10)  fixed,  are: 

uh(x,z)  =  6S 1 1  (L2-x2)z+2S  ,  3z3+S55(2z3-  |h2z)1 

bH3L  2  (A8) 


wh(x,z)  =  — — ^2Sn(x3-3L2x+2L3)  -  6Sj3Xz2J 

where,  L,b,  and  H  are  the  length,  width,  and  depth  of  both  beams,  respectively.  The  effective 
compliance  matrix  for  the  homogeneous  beam,  Sjj  (i,j=l,3,5),  is  to  be  determined  by  matching 

the  boundary  displacement  of  the  homogeneous  beam  with  that  of  the  laminated  beam. 
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In  Problem  1,  the  elements  of  the  effective  compliance  matrix,  Sy,  were  predicted  by  matching 
the  outer  and  inner  displacements  of  the  laminated  and  homogenized  rings  exactly.  However  in 
this  problem  the  boundary  displacement  field  of  the  homogeneous  beam  is  matched  with  that  of 
the  laminated  beam  at  a  discrete  number  of  boundary  points  by  minimizing  the  error.  A: 


n 

A  =  X  [  {  Uh(Xi,Zi)-Uc(Xi,Zi)  }  2+  {  Wh(Xi,Zi)- Wc(Xi,Zi)  }  2] 
i=l 


(A10) 


where  \[,z[  are  coordinates  of  the  points  selected  along  the  boundary  of  the  two  beams. 


After  having  estimated  the  elements  of  Sy  by  minimizing  A  in  equation  (A  10),  the  effective 
bending  stiffness,  and  interlaminar  shear  stiffness  and  Poisson's  ratio  of  an  orthotropic  and 
symmetric  laminate  for  bending  mode  of  deformation  are  respectively  : 


(All) 


The  stress  analysis  in  a  similar  fashion  for  stresses  in  the  2-3  plane  yields  the  following 
additional  effective  elastic  constants: 


1 


~f  -  S23 

V32"  S' 


G  23  ~  c  „  a 

b44  ^  22  (A  12) 

A  sensitivity  study  was  performed  on  the  convergence  of  the  values  of  the  effective  properties 
by  varying  the  number  of  points  on  the  boundary.  It  was  found  that  the  bending  stiffness  and 
the  interlaminar  shear  stiffness  converged  very  fast  with  increasing  numbers  of  points  on  the 
boundary  and  were  practically  insensitive  to  the  ratio  of  the  number  of  points  chosen  on  the 
vertical  lines  to  that  of  the  horizontal  lines  of  the  boundary.  The  interlaminar  Poisson’s  ratios 
were,  however,  somewhat  sensitive  to  the  latter  ratio.  We  found  out  that  the  value  of  the  ratio 
equal  to  one  gave  a  good  estimate  of  the  Poisson’s  ratio,  thus  this  value  of  the  ratio  was  used  in 
obtaining  the  results  presented  later. 


Problem  3:  Beam  with  Uniformly  Distributed  Load 

The  beam  considered  is  under  uniformly  distributed  load  acting  on  the  top  and  bottom  surfaces 
of  the  beam  of  intensity  q  and  p,  respectively,  (Figure  1 1).  Besides  all  the  effective  properties 
obtained  from  Problem  2,  one  more  additional  property,  the  interlaminar  effective  normal 
rif 

stiffness,  ^3,  is  obtained  from  this  problem.  The  method  of  predicting  the  effective  moduli  in 
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this  problem  is  similar  to  that  of  Problem  3  and  thus  is  not  discussed.  For  details  one  may  refer 
to  reference  [7]. 
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Figure  1 1 .  Configuration  of  the  laminated  beam  and  its  equivalent  homogeneous  beam  under 
uniformly  distributed  load  on  the  top  and  bottom  surface. 


In  general,  the  effective  moduli  of  a  laminate  depend  on  the  laminate  stacking  sequence  and  mode 
of  deformation.  When  the  laminate  contains  a  large  number  of  sublaminates,  the  effect  of 
stacking  sequence  almost  disappears;  then  the  effective  moduli  become  practically  independent  of 
mode  of  deformation.  The  values  of  effective  moduli  given  below  are  independent  of  mode  of 
deformation,  thus  the  superscripts  'o'  or  'f  are  omitted  from  the  notations  below. 


1 )  Three-dimensional  properties  of  T300/5208  composites:  [0]  (unidirectional  lamina) 
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2)  Three-dimensional  effective  properties  of  T300/5208  quasi-isotropic  laminate,  |tt/4]s 
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0. 

0. 

1 

P 

E,  u»t2 

u13 

69.7 

0.30 

0.03 

C55 

0 

i 

5.7 

0. 

u2\  E2 

y23 

= 

0.30 

69.7 

0.03 

C66 

GPa 

26.8 

U31  l/3  2 

E3 

0.15 

0.15 

12.9 

C44 

=  S44-1  = 

G23. 

C55  =  S55 

1  =  G31-  C66  =  S66 

1  =e12 

3)  Three-dimensional  properties  of  IM6/Epoxy  composites:  [0]  (unidirectional  lamina) 


r  r  r 

°xx  uxz 

2089  9.2  9.2 

s  s  s 

°xx  °xy  Jv 

4.9  -1.6  -1.6 

c  c 

°yy  yz 

= 

17.9  10.9 

Gyy  Syz 

= 

89.3  -53.6 

1 

0 

N 

N 

l _ 

GPa  17.9 

SrZ 

TPa-’  89.3 

O 

jO 

O 

O 

_ l 

3.5 

0. 

0. 

-  -1 

h- 

I'xy 

^xz 

203.0  0.02 

0.02 

crr  0 

= 

8.4 

0. 

^yx 

V 

= 

0.32 

1  1.2 

0.60 

CSS 

GPa 

8.4 

/zx 

^zy 

Ez 

0.32 

0.60 

1  1.2 

Cqq  = 

5  -1 

Oqq 

II 

m 

r  -  s 

^rr  “  °rr 

1  =  E 

r>  Gss 

-  S 

-  °ss 

=  ES 
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4)  Three-dimensional  effective  properties  of  IM6/Epoxy  quasi-isotropic  laminate,  [tc/4]s 


Eli  C12 

E13 

87.1 

26.7 

4.41 

<N 

co" 

in 

Sis' 

12.8 

-3.8 

-2.6 

C22 

C23 

= 

87.1 

4.4 

S22 

S23 

- 

12.8 

-2.6 

^33 

GPa 

15.4 

S33 

TPa"1 

66.3 

+ 


C44  0  0 

5.0 

0. 

0. 

"e, 

u12 

u13 

78.4 

0.30 

0.04 

S55  0 

= 

5.0 

0. 

y2, 

E2 

u23 

= 

0.30 

78.4 

0.04 

C66 

GPa 

30.2 

1^31 

^3  2 

E3 

0.20 

0.20 

15.1 

C44  -  S44  ’  =  G23 ,  C55  =  S55  1  =  G31,  C66  =  See"1  =  G12 


5)  Three-dimensional  properties  of  E-Glass/Epoxy  composites:  [0]  (unidirectional  lamina) 


r  r  r 

uxx  ^xy  uxz 

41.6  5.8  5.8 

Sxx  Sxy  5XZ 

25.9  -6.7  -6.7 

C  C 

*~yy  uyz 

zz 

13.7  8.6 

s 

'-’yy  wyz 

= 

120.9  -72.6 

c“. 

GPa  13.7 

. 

S“. 

TPa-’  120.9 

+ 


cqq  0  0  ‘ 

2.6 

0. 

0, 

"Ex 

I'xy 

yxr 

38.6  0.06 

0.06 

Crr  0 

= 

4.2 

0. 

>V 

Ey 

Vvz 

= 

0.26 

8.3 

0.60 

Css 

GPa 

42 

N 

Ez 

0.26 

0.60 

8.3 

f  -S' 
-  °qq 


=  E, 


r  -  s 

^rr  -  Jrr 


-1 


Er/  Ess  =  2Ss  1  =  Es 


6)  Three-dimensional  effective  properties  of  E-Glass/Epoxy  quasi-isotropic  laminate, 
lrc/4  is 
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Cl,  c12  C,3 

22.3  7.4  4.8 

Si  1  S12  s13 

52.7  -14.2  -15.0 

C22  C23 

- 

22.3  4.8 

S22  S23 

= 

39.9  -15.0 

C33 

GPa  12.3 

S33 

TPa'1  93.3 

+ 


- 1 

n 

a 

A 

O 

O 

_ l 

3.2 

0. 

0. 

E, 

y12 

1^13 

1 9.0 

0.27 

0.16 

C55  0 

= 

3.2 

0. 

U2i 

C2 

v23 

= 

0.27 

19.0 

0.16 

Ce6 

GPa 

7.5 

U31 

u3  2 

E3 

0.29 

0.29 

10.7 

C44  =  S44  1  =  G23 ,  c55  =  s55  1  =  g31,  c66  =  S66  1  -  g12 


)\ 
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Appendix  B 


INSTRUCTIONS  FOR  EXECUTING  THE  FORTRAN  PROGRAM  "sphere" 

To  execute  the  code  one  needs  to  open  the  application  file  "sphere  apl",  then  sequentially  follow  the 
input-output  sequence.  The  code  uses  a  material  data  file  "matdat.spr".  This  file  contains  the 
effective  stiffness  properties  of  quasi-isotropic  laminates,  and  strength  properties  and  Young’s 
moduli  of  the  unidirectional  plies.  The  data  format  of  the  file  is  printed  on  a  following  page.  One 
can  easily  edit  the  file  to  modify  the  built-in  data  or  add  new  data.  The  method  of  calculating  the 
three-dimensional  effective  stiffnesses  is  presented  earlier  in  Appendix  A.  A  printout  of  the  input- 
output  sequence  of  an  execution  run  of  the  code  is  also  included  in  this  report. 


sphere.for 


Thu,  Sap  6, 1990 
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PROGRAM  SPHERE 

r* 

c —  This  procedure  is  to  do  stress  and  failure  analyses  of  a  thick 

c —  multilayered  sphere  subjected  to  INTERNAL,  EXTERNAL  pressures, 

c —  Each  layer  is  assumed  to  be  composed  of  quasi-isotropic  lay  ups, 

c —  each  quasi-isotropic  layer  is  thus  assumed  to  behave  as 

c —  transversely  isotropic  layer  (isotropic  plane  is  theta-phi 

c —  plane,  perpendicular  to  the  radial  direction) .  The  code  reads 

c —  Material  data  from  an  existing  data  file  named  'matdat . spr ' . 

c —  The  material  data  are  given  in  'SI'  units.  The  material 

c —  data  represent  the  three-dimensional  effective  properties  n* 

c —  a  quasi-homogeneous  lay  up.  The  format  for  the  material  data 

c —  file  is  given  below.  The  field  length  for  numenric  data:E9.3 


c — 

number  Name 

Cll 

C12 

C13 

C22 

C23 

c — 

C33 

C44 

C55 

C66 

X 

c — 

X' 

Y 

Y' 

S 

Fxy 

c — 

Fyz* 

alpl 

alp2  bta2 

Em* 

c — 

El 

E2 

nu21  nu32 

Es 

c —  The  geometric  and  lay-up  configurations  can  be  given  interactively 

c —  which  should  be  consistent  with  the  unit  used  in  the  data  file, 

c —  'matdat . spr ' .  A  quadratic  failure  criterion  (TSAI-WU)  is  used  for 

c —  predicting  failure  loads.  The  procedure  is  limited  to 

c —  analyzing  rings  (pressure  vessels)  of  maximum  25  layers. 

c —  developed  by 


c — 

★ 

A JIT  K.  ROY 

★ 

c — 

★ 

UDRI 

★ 

c — 

★ 

300  College  Park 

★ 

c — 

★ 

DAYTON,  OH  45469-0001 

★ 

c — 

★ 

(513)  255-9104 

★ 

C 


REAL  KM 
CHARACTER* 9 
CHARACTER* 1 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 


DF1 

TAB 

SM11 (25) , SMI 2 (25) , SM13 (25) , SM33 (25) ,  KM  (25) , CM (25) 
QUE(26) ,SGR (25, 11) ,SGT (25,11 ) ,EU (25,11 ) ,CT11 (25) 
CT12 (25) , CT13 (25) , CT33 (25) , AR (25, 11 ) , PROP (10,25) 
STRR (25,11) , STRT (25,11) 

AM (26) ,YY (25,11 ) ,INDX (25) 


EXTERNAL  F, FPRIM 


COMMON  / CLK/  SM11, SM12, SM13, SM33, KM, CM, AM 
COMMON  / STRNTH/  CT11, CT12, CT13, CT33, PROP 
COMMON  / SSDP/  SGR, SGT, STRR, STRT, EU,YY,AR 


c —  the  subroutine  'PLYLUP'  is  for  sphere  geometry,  material 
c  property  input,  and  calculations  for  layer  stiffness  and 
c  compliance. 


CALL  PLYLUP (INDX, QUE, NM, DF) 


c—  the  subroutine  'CALQUS'  is  for  the  calculation  of  the 
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c  interlayer  tractions  'q* 

CALL  CALQUS (NM, QUE ) 

c —  the  subroutine  'STSDSP'  is  for  the  computation  r  the 

c  stresses,  strains,  and  radial  displacement 

CALL  STSDSP (NM, QUE,  ND) 


WRITE ( * , * )  'INPUT  OUTPUT  FILE  NAME  (max  9  char)' 
READ (*, ' (A) ' )  DF1 
OPEN ( 7, FILE=DF1,STATUS=' UNKNOWN ' ) 

WRITE (*,*)  '***  PLEASE  WAIT  ***' 


TAB  =  CHAR (9) 

C 

WRITE (7,*)  'inner  radius (meter) ,  outer  radius (meter ) ' 

WRITE (7,*)  AM ( 1 ) , TAB , AM (NM+ 1 ) 

WRITE (7,*)  'b/a  SPHERE' 

BOA  =  AM ( NM+ 1 ) /AM ( 1 ) 

WRITE (7,*)  BOA 

WRITE (7,*)  'mat  #  r, meter  sig(r),MPa  sig (th) =sig (phi) , 

+MPa' 

WRITE (7, 340)  ( (INDX (I ) , TAB, YY (I, J) , TAB, 

+  SGR(I, J) /l.E6,TAB,SGT(I, J)/1.E6, J=1,ND) ,I=1,NM) 
write (7,*)  'mat  #  r, meter  eps(r)  eps (th)=eps (phi) 

+u (r) ,mm' 

WRITE (7, 345)  ( (INDX (I ) , TAB, YY (I, J) , TAB, 

+  STRR (I, J) , TAB, STRT (I, J) , TAB, EU (I, J) *1000, J=1,ND) , 1=1, NM) 

C 

C 

CALL  STRIO (NM, ND, ITK, JTK, ARMN, DF, INDX) 

C 

WRITE (7,*)  'mat  #  r (meter)  strength  ratio' 

WRITE (7, 342)  ( (INDX (I) , TAB, YY (I, J) , TAB, AR (I, J) , J=l, ND) , 1=1, NM) 

WRITE (7,*)  'FAILURE  LOCATION  &  FAILURE  PRESSURE' 

WRITE (7, *)  'mat  #  r(meter),  INT  PRESSURE,  EXT  PRESSURE  (in  Pa)' 
WRITE (7, 345)  INDX (ITK) , TAB, YY (ITK, JTK) , TAB, AR (ITK, JTK) *QUE(1) , 

+  TAB, AR (ITK, JTK) *QUE (NM+1 ) 

WRITE (7, 360) 

360  FORMAT (’END  OF  OUTPUT') 

340  FORMAT (2X,I2,A1,E12.6,A1,E12.6,A1,E12.6) 

342  FORMAT (2X,I2,A1,E12.6,A1,E12.6) 

345  FORMAT (2X,I2,A1,E12.6,A1,E12.6,A1,E12.6,A1,E12.6) 

CLOSE ( 7 , STATUS= ' KEEP ' ) 

STOP 

END 

r- 

r  ***★***•*********★★*★**★**★*★*★★**★★*****★★*★★★**★ 

SUBROUTINE  PLYLUP (INDX, QUE, NM, DF) 


This  subroutine  is  to  obtain  the  geometric  and  lay-up 
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c —  configuration  of  the  ring.  The  Ply  Degradation  Factor  (DF) 

c —  is  also  given  interactively  for  the  prediction  of  Last-Ply-Fail 

c —  -ure  or  Burst  Pressure. 

c 

CHARACTER  TYP*21,  WISH*7 
REAL  NPLY, KM 

DIMENSION  PROP (10,25) , AM (26) , PLYPRP (25) , SM11 (25) ,SM12 (25) 
DIMENSION  SM13 (25) ,SM33(25) ,  KM  (25) ,  CM(25) 

DIMENSION  CTl 1 (25) ,CT12(25) ,CT13(25) ,CT33(25) 

DIMENSION  INDX(25) , NPLY (25)  ,  QUE (26) 

c 

COMMON  / CLK/  SM11 , SM12, SMI 3, SM33, KM, CM, AM 
COMMON  /STRNTH/  CT11 , CT12, CTl 3, CT3 3, PROP 


WRITE (*,21) 

21  FORMAT (/, ’INPUT  INNER  DIAMETER',/, '?') 

READ ( * , * )  DMT I NR 

C  DMTINR=2 . 0 

C  LNTH=1 . 0 

C 

WRITE (*,*)  'INPUT  THE  VALUE  OF  b/a  OF  THE  SPHERE' 

READ ( * , * )  BOAS 

C 

WRITE (*,22) 

22  FORMAT (’ INPUT  INTERNAL  PRESSURE  &  EXTERNAL  PRESSURE',/,’?') 
READ ( * , * )  QUE1,QUEN 

C  QUE 1-1 . 0 

C  QUEN= i . 0 

C 
C 

OPEN ( 8 , F I LE= ' mat dat . spr ' , STATUS- ' OLD ' ) 

C 

WRITE (*, 10) 

10  FORMAT (/, T25, ' ****MATERIAL  INFORMATION****') 

WRITE (*, 11) 

11  FORMAT (/,T1 5, 'INDEX  NUMBER  MATERIAL  DESCRIPTION' , 

+  /, T15,  ' - ') 

C 

1  =  0 

1  CONTINUE 
I  =  1  +  1 

READ  (8,20,  END=2 )  IND,  TYP,  (PROP  (I,  J) ,  J=l,  25) 

20  FORMAT (2X,I2,A21,T27,5E9.3,/,T27,5E9.3,/,T27,5E9.3,/,T27,5E9.3 

+  , /,T27, 5E9. 3) 

C 

WRITE (*, 12)  IND, TYP 
1?  FORMAT (T23, 12,T30, A2  3 ) 

GOTO  1 

C 

2  1=1-1 

CLOSE (UNIT-8,  STATUS- ' KEEP ' ) 

ITR  =  0 

C  WRITE (*,*)  'input  MATERIAL  NUMBER  for  the  mandrel’ 
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C  READ ( * , * )  INDX(l) 

WRITE (*,*)  '***  INPUT  PLY  LAYUP  INFORMATION  (max  25  layers) 
C51  IR=IR+1 

WRITE (*,*) 

C 

C  WRITE (*,*)  'INPUT  THE  VALUE  OF  THE  DEGRADATION  FACTOR' 

C  WRITE (*,*)  ’  for  INTACT  plies:  DF=1' 

C  READ ( * , * )  DF 

DF=1 . 0 

C  BOAS=  l.+l./AOTS 


c  mandrel  displacement 

AM(1)  =  0. 5*DMTINR 

C  AM (2)  =  BOAM*AM ( 1 ) 

C  NPLY(l)  =  AM(2)  -  AM ( 1 ) 

WRITE (*,*)  'INPUT  NO  OF  LAYERS  (max  25  layers)’ 
READ ( * , * )  NM 


DO  25  1=1 , NM+1 
QUE(I)  =  0.0 
25  CONTINUE 

QUE(l)  =  QUE1 
QUE (NM+1 )  =  QUEN 

DO  51  IR=1 , NM 

WRITE (*,*)  '**  LAYERS  ARE  NUMBERED  FROM  INSIDE  TO  OUTSIDE  * 
WRITE (*,*) 

WRITE (* , * )  '***  FOR  LAYER  ' , IR, '  ***’ 

WRITE (*,*) 'INPUT  MATL  INDEX  NO.  &  THICKNESS  FACTOR' 

C 

READ (*,  *)  INDX(IR)  ,NPLY(IR) 

NPLY (IR) =NPLY (IR) * (BOAS-1. ) *AM(1) 

51  CONTINUE 


DO  53  IR=1 , NM 

DO  50  I  -  1,25 

PLYPRP ( I )  =  PROP (INDX (IR) ,  I) 

50  CONTINUE 

C 

C —  degradation  of  effective  stiffness  and  compressive  strength 

PLYPRP (2)  =  DF*  PLYPRP (2) 

PLYPRP ( S )  =  DF*PLYPRP(3) 

PLYPRP (5)  =  DF*  PLYPRP (5) 

C 

AM ( IR+1 )  =  AM ( IR)  +  NPLY ( IR) 

CM ( IR)  =  AM (IR) /AM ( IR+1 ) 

C 

CALL  BTCLCT (PLYPRP, IR) 
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C 

KM(IR)  =  0 . 5*SQRT (1 . +8 . * (CT11 (IR) +CT12 (IR) -CT13 (IR) ) /CT33 (IR) ) 

WRITE (*,*)  'DO  YOU  WANT  TO  ADD  ANOTHER  PLY  LAYUP?  (Y  OR  N) ' 
READ (*, ' (A) ' )  WISH 


3  CONTINUE 

IF (WISH  .EQ.  'Y')  GOTO  51 


RETURN 

END 


*************************************************** 
SUBROUTINE  BTCLCT (PLYPRP, ITR) 


c —  this  subroutine  is  to  perform  some  intermediate  calculations 

c —  required  for  the  routine  'CALQUS'. 

c 

REAL  KM 


c 


C 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


PLYPRP (25) , SMTX (6, 6) , CT11 (25) , CT12 (25)  ,  CT13 (25) 
CT33 (25) , AM (26) , CTX (6, 6) ,  PROP (10,25) 
XX(6,1),KM(25),CM(25) 

SM11 (25) , SMI 2 (25) , SM13 (25) , SM33 (25) 


COMMON  /CLK/  SM11 , SM12, SMI 3, SM33, KM, CM, AM 
COMMON  /STRNTH/  CT11 , CT12, CT13, CT33, PROP 


11 

10 

c 

c 

c 


DO  19  I  —  1,6 
XX(I,1)  =  0.0 
DO  11  J  =  1,6 
CTX (I, J)  =  0.0 
SMTX (I, J)  =0.0 
CONTINUE 
CONTINUE 


CTX (1,1) 
CTX (1,2) 
CTX  (1,3) 
CTX (2,1) 
CTX (2, 2) 
CTX (2, 3) 
CTX  (3,1) 
CTX (3, 2) 
CTX ( 3, 3) 
CTX (4,4) 
CTX (5, 5) 
CTX (6, 6) 


PLYPRP (1) 
PLYPRP (2) 
PLYPRP (3) 
CTX (1,2) 
PLYPRP (4) 
PLYPRP (5) 
CTX  (1,3) 
CTX (2, 3) 
PLYPRP (6) 
PLYPRP (7) 
PLYPRP (8) 
PLYPRP (9) 


DO  15  1=1,6 
DO  16  J=l,  6 
SMTX ( I , J)  =  CTX (I, J) 
16  CONTINUE 
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15  CONTINUE 

CALL  GAUSS J (SMTX, 6,6, XX, 1,1) 

C 

c 

CTll(ITR)  =  CTX  (1,1) 

CT12 (ITR)  =  CTX (1,2) 

CT13 (ITR)  =  CTX (1,3) 

CT33 (ITR)  =  CTX (3, 3) 
c 

SM11 (ITR)  =  SMTX (1, 1) 

SMI 2 (ITR)  =  SMTX (1,2) 

SMI  3 (ITR)  =  SMTX (1,3) 

SM33 (ITR)  =  SMTX (3, 3) 

C 

RETURN 

END 

c  *★★★★★★*★*★**★★*★*★★**★*★★★★**★★★***★★★***★★****★★*•★ 

SUBROUTINE  CALQUS (NM, QUE) 

c —  this  subroutine  is  to  calculate  the  inter-layer  normal  pressure 
c —  Q's  from  a  set  of  simultaneous  equations. 

REAL  KM,  MUE 
REAL  AA, BB 

DIMENSION  SM11 (25)  ,SM12 (25) , SM13 (25) , SM33 (25) ,  CM (25) 

DIMENSION  AA (25, 25) , QUE (26) , ALF (25) , BTM(25) , BB (25, 1) 

DIMENSION  SGR (25, 11) , SGT (25,11) , STRR (25,11) ,  PROP (10,25) 
DIMENSION  STRT (25, 11 ) , EU (25,11), YY (25,11), AR (25, 11 ) 

DIMENSION  CT11 (25) ,CT12 (25) , CT13 (25) ,CT33 (25) 

DIMENSION  AM (26) , KM (25) ,MUE (25) , GMA (25) 

COMMON  /CLK/  SM11, SM12, SM13, SM33, KM, CM, AM 
COMMON  / STRNTH/  CT11, CT12, CT13, CT33, PROP 
COMMON  /SIZE/  AA, BB 

COMMON  / SSDP/  SGR, SGT, STRR, STRT, EU,YY,AR 


DO  10  I  =  1,  NM-1 

BB (I,  1)  =  0.0 

DO  11  J  =  1 , NM-1 
AA (I, J)  =  0.0 
11  CONTINUE 

10  CONTINUE 

DO  20  1=1, NM 

MUE (I )  =  (1 . -CM (I ) ** (2 . *KM (I ) ) ) * ( (KM (I ) + . 5) *CT33 (I ) 

1  -2.*CT13(I) )*( (KM(I)-.5) *CT33 (I) +2 . *CT13 (I) ) 

ALF ( I )  =  2 . *KM (I ) *CT33 (I) *CM (I ) ** (KM (I) + . 5) /MUE (I) 


BTM(I)  =  ( 1 . -CM (I ) ** (2 . *KM (!)))*(. 5*CT33 (I ) -2 . *CT13 (I ) ) 
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GMA(I)  =  KM  (I )  *  (1 .  +CM  (I )  **  (2 .  *KM  (I )  )  )  *CT33  (I )  /MUE  (I ) 

20  CONTINUE 

IF  (NM  .EQ.  1)  GOTO  50 
IF ( (NM-1 )  .EQ.  1)  THEN 

AA(1, 1)  =  (BTM(2) -BTM(l) -GMA (2) -GMA (1) ) *AM (2) 

BB ( 1 , 1 )  =  -QUE (NM+1 ) *AM (3) *ALF (2) /CM (2) -QUE (1 ) *AM (1 ) *ALF (1 ) 
QUE (2)  =  BB (1) /AA(1, 1) 

ELSE 

DO  40  1=1, NM-1 

IF (I  .EQ.  1)  THEN 

BB  ( 1 , 1 )  =  -AM ( I ) *ALF ( I ) *QUE ( 1 ) 

AA  (1,1)  =  AM(2)*(BTM(I+1)-BTM(I)-GMA(I+1)-GMA(D) 

AA (I, 1  +  1 )  =  AM (3) *ALF (I-H  ) /CM(I+1) 

ELSEIF  (I  .EQ.  (NM  THEN 

BB  (1,1)  =  -AM  (H  *ALF  (1+1)  *QUE  (NM+1)  /CM(I+1) 

AA  (I ,  NM-1 )  =  AM(I+1)*(BTM(I+1)-BTM(I)-GMA(I+1)-GMA(I)  ) 

AA (I,NM-2)  -  AM(I) *ALF (I) 

ELSE 

BB (1,1)  =  0.0 

AA(I, 1-1)  =  AM (I ) *ALF (I ) 

AA(I, I)  =  AM (1+1) * (BTM (1+1) -BTM (I) -GMA (1+1) -GMA (I) ) 

AA (I, 1+1 )  =  AM(I+2)*ALF(I+1) /CM(I+1) 

ENDIF 

40  CONTINUE 

NMM1  =  NM-1 
IFAIL  =  0 

CALL  GAUSSJ (AA, NMM1, 25,  BB,  1,1) 

DO  30  1=1,  NM-1 
QUE ( I + 1 )  =  BB (1,1) 

30  CONTINUE 

ENDIF 

50  CONTINUE 


RETURN 

END 

Q  ***************************************************** 

SUBROUTINE  STSDSP (NM, QUE, ND) 
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c —  to  calculate  stress,  strain  and  displacement  field  based  on 

c —  the  thick  sphere  formulation. 

REAL  KM 

DIMENSION  YY(25,11) , SGR (25,11 ) , SGT (25,11 ) ,EU (25,11 ) ,SM11 (25) 
DIMENSION  SMI 2 (25) ,SM13 (25) ,SM33 (25) , CM (25) , AM(26) , KM(25) 
DIMENSION  CT11 (25) , CT12 (25) , CT13 (25) , CT33 (25)  ,  PROP (10, 25) 
DIMENSION  AR(25, 11) , STRR(25, 11)  ,STRT(25,11) 

DIMENSION  QUE (26) 

COMMON  / CLK/  SMll , SM12, SM13, SM33, KM, CM, AM 
COMMON  / STRNTH/  CT11, CT12, CT13, CT33, PROP 
COMMON  /SSDP/  SGR, SGT, STRR, STRT,  EU,  YY,  AR 

WRITE (*, *) 

WRITE (*,*)  'INPUT  NO  OF  EQUIDISTANT  CALCULATION  POINTS  THRU 
+THICKNESS  IN  EACH  LAYER' 

WRITE (*,*)  '(max  11,  min  2)' 

READ  ( * , * )  ND 

DO  40  I  =1,NM 

FCT  =  (AM (1+1) -AM (I) ) / (ND-1) 

CNST1  =  1./ (l.-CM(I)** (2.*KM(I) ) ) 

CNST2  =  (KM (I) - . 5) *CT33 (I)  +  2.*CT13(I) 

CNST3  =  (KM (I) + .  5) *CT33 (I)  -  2.*CT13(I) 

CNST4  =  CT11 (I ) +CT12 (I ) + (KM(I ) - . 5) *CT13 (I) 

CNST5  =  CT11 (I) +CT12 (I) - (KM (I) + . 5) *CT13 (I) 

DO  45  J  =  1 , ND 

YY (I,  J)  =  AM(I)  +  ( J-l ) *FCT 

c —  I  :  for  layers,  increasing  from  inside  to  outside 

c —  J  :  for  calculation  points,  incresing  from  inside  to  outside 


c —  radial  stress 

SGR (I, J)  =  CNSTl* ( (QUE (I) *CM(I) ** (KM (I) +1.5) -QUE (1+1) ) * 

1  (YY(I,J)/AM(I+1) )**(KM(I)-1.5)  - 

2  (QUE (I) -QUE (1+1) *CM(I) ** (KM (I) -1.5) )* 

3  CM (I) ** (KM (I) +1.5) * (AM (I+1)/YY(I,J))** (KM (I) +1.5) ) 

c —  hoop  stress 

SGT (I, J)  =  CNSTl* (CNST4/CNST2) * (QUE (I) *CM(I) ** (KM(I) +1 . 5) 

-QUE (1+1) ) * (YY (I, J) /AM(I+1) ) ** (KM (I) -1.5)  + 

CNSTl* (CNST5/CNST3) * (QUE (I) -QUE (1+1) * 

CM  (I) ** (KM(I)-l . 5) ) * 

4  CM (I)** (KM (I) +1.5)* (AM(I+1) /YY(I, J) ) ** (KM (I ) +1 . 5) 


45  CONTINUE 

40  CONTINUE 


DO  50  1=1, NM 
DO  60  J=1 ,  ND 
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c —  radial  strain 

STRR (I, J)  =  SM33 (I) *SGR (I, J) +2 . *SM13 (I ) *SGT (I, J) 
c —  hoop  strain 

STRT (I, J)  =  SMI 3 (I ) *SGR(I, J)  +  (SM11 (I ) +SM12 (I ) ) *SGT (I,  J) 

c —  radial  displacement 

EU (I, J)  =  STRT (I,  J) *YY (I, J) 

60  CONTINUE 
50  CONTINUE 

RETURN 

END 

SUBROUTINE  STRIO (NM, ND, ITK, JTK, ARMN, DF, INDX) 

c —  to  obtain  the  strength  ratio  for  predicting  the  failure 

c —  pressure. 

DIMENSION  SM11 (25) , SM12 (25) , SM13 (25) , SM33 (25) , KM (25) 
DIMENSION  CM (25) ,AM(25) , INDX(25) , PROP (10, 25) , CMTX (6, 6) 
DIMENSION  CT11 (25) , CT12 (25) , CT13 (25) , CT33 (25) 

DIMENSION  SGR(25, 11) , SGT (25, 11) , STRR (25, 11) , STRT (25,11) 
DIMENSION  EU (25,11), YY (25,11) , AR (25, 11 ) , PLYPRP (25) 
DIMENSION  GE (6, 6) , GEE (6) , AF (6, 6) ,  AFF (6) 

COMMON  /CLK/  SMll , SM12, SM13, SM33, KM, CM, AM 
COMMON  /STRNTH/  CT11, CT12, CT13, CT33, PROP 
COMMON  /SSDP/  SGR, SGT, STRR, STRT,  EU,YY,AR 

DO  2  1=1,6 

AFF  (I)  =  0.0 

DO  3  J=l, 6 

AF (I, J)  =  0.0 

3  CONTINUE 

2  CONTINUE 


DO  10  1=1, NM 

DO  50  IR  =  1,25 

PLYPRP (IR)  =  PROP (INDX (I) , IR) 

50  CONTINUE 

c —  degradation  of  unidirectional  stiffness  properties 

PLYPRP (11 )  =  DF**0. 2*PLYPRP (11) 

PLYPRP (22)  =  DF*  PLYPRP (22) 

PLYPRP (23)  =  DF*PLYPRP (23) 

PLYPRP (25)  =  DF*  PLYPRP (25) 


c —  3:  fiber  direction 
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c —  1  &  2:  transverse  to  the  fiber 

AF (1,1)  =  1./ (PLYPRP (12) *PLYPRP (13)) 

AF (2,2)  =  AF (1,1) 

AF (3, 3)  =  l./(PLYPRP(10)*PLYPRP(ll)) 

AF (4,4)  =  1./ (PLYPRP (14 ) **2) 

AF  (5, 5)  =  AF (4,4) 

AF  (6, 6)  =  0. 60/PLYPRP (10) 

AFF(l)  =  1. /PLYPRP (12)  -  1 . /PLYPRP (13) 

AFF  (2)  =  AFF (1) 

AFF ( 3)  =  1. /PLYPRP (10)  -  1. /PLYPRP (11) 

AF (1,2)  =  PLYPRP(16)*SQRT(AF(1,1)*AF<2,2)) 
AF (1,3)  =  PLYPRP (15) *SQRT(AF (1, 1) *AF (3, 3) ) 
AF  (2,3)  =  AF  (1,3) 

AF (2,1)  =  AF (1,2) 

AF (3,1)  =  AF (1,3) 

AF  (3,  2)  =  AF  (2,  3) 


CALL  CTX 3D (PLYPRP, CMTX) 

CALL  G3DCLT(CMTX,AF,AFF,GE,GEE) 
DO  11  J=1 , ND 


AA  =  ABS ( (GE (3, 3) +2 . *GE (3,2) +GE (2,2) ) *STRT (I, J) **2 

1  +2 . * (GE (3, 1 ) +GE (2, 1 ) ) *STRT (I,  J) *STRR (I,  J) 

2  +GE (1,1) *STRR(I, J) **2) 

BB  =  (GEE (3) +GEE (2) ) *STRT (I, J) +GEE (1 ) *STRR (I,  J) 

c —  strength  ratios  based  on  thick  preeuse  vessel  solution 

AR  (I , J)  =  - (BB/ (2 . *AA) )  +SQRT ( (BB/ (2 . *AA) ) **2  +1./AA) 


11  CONTINUE 

10  CONTINUE 

c —  ARMN :  minimum  value  of  strength  ratios  based  on  thick  pressure 

c —  vessel  solution. 


CALL  VARMN (AR, ARMN, NM,  ND, ITK,  JTK) 

RETURN 

END 

★  ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★■AT*-*"***'* 

SUBROUTINE  VARMN (AR, ARMN, NM, ND, ITK, JTK) 

DIMENSION  AR (25, 11 ) 

ITK  =  1 
JTK  =  1 

ARMN  =  AR (1,1) 
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DO  12  1=1, NM 
DO  11  J  =  1 , ND 

IF (AR(I, J)  .GE.  ARMN)  GOTO  11 

ITK  =  I 

JTK  =  J 

ARMN  =  AR(I, J) 

11  CONTINUE 

12  CONTINUE 

RETURN 

END 

SUBROUTINE  G3DCLT (CTX, AF, AFF, GE,  GEE) 

DIMENSION  CTX (6, 6) , GE (6, 6) ,TMP1 (6,6), GEE (6) , TMP2 (6) 
DIMENSION  AF (6,6),  AFF (6) 

DO  5  1=1,6 
DO  7  J=l, 6 
TMP2(I)  =0.0 
TMP1 (I, J)  =  0.0 
GE (I, J)  =0.0 
7  CONTINUE 

5  CONTINUE 

DO  10  1=1,6 
DO  11  J=l, 6 
DO  12  K=l, 6 

TMP1  (I,  J)  =  TMP1  (I ,  J)  +  AF(I,K)*CTX(K,  J) 

12  CONTINUE 

11  CONTINUE 

10  CONTINUE 

DO  20  1=1,6 
DO  21  J=l, 6 
DO  22  K=l, 6 

GE (I, J)  =  GE (I, J)  +  CTX (I, K) *TMP1 (K, J) 

22  CONTINUE 

21  CONTINUE 

20  CONTINUE 

DO  30  1=1,6 
DO  31  J=l, 6 

TMP2 (I)  =  TMP  2(1)  +  AFF ( J) *CTX ( J, I ) 

31  CONTINUE 

30  CONTINUE 

DO  40  1=1,6 
GEE (I)  =  TMP2 (I ) 

40  CONTINUE 

RETURN 

END 

c —  the  routine  for  solving  a  system  of  simultaneous  equations 
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SUBROUTINE  GAUSSJ(A,N,NP,B,M,MP) 

PARAMETER  (NMAX=50) 

DIMENSION  A (NP, NP) , B (NP, MP) , IPIV (NMAX) , INDXR (NMAX) , INDXC (NMAX) 
DO  11  J=1,N 
IPIV(J)=0 

11  CONTINUE 

DO  22  1=1, N 
BIG=0. 

DO  13  J=l, N 

IF (IPIV ( J) . NE . 1) THEN 
DO  12  K=1,N 

IF  (IPIV (K) . EQ. 0)  THEN 

IF  (ABS (A ( J,  K) )  . GE . BIG) THEN 
BIG=ABS (A ( J, K) ) 

IROW=J 

ICOL=K 

ENDIF 

ELSE  IF  (IPIV (K) . GT . 1 )  THEN 
PAUSE  'Singular  matrix* 

ENDIF 

12  CONTINUE 
ENDIF 

1 3  CONTINUE 

IPIV (ICOL) =IPIV (ICOL) +1 
IF  (IROW.NE. ICOL)  THEN 
DO  14  L=-1,N 
DUM=A ( IROW, L ) 

A (IROW, L) =A (ICOL, L) 

A (ICOL, L) =DUM 

14  CONTINUE 

DO  15  L=1,M 
DUM=B (IROW, L) 

B (IROW, L)=B (ICOL, L) 

B (ICOL, L) =DUM 

15  CONTINUE 
ENDIF 

INDXR (I ) =IROW 
INDXC ( I ) =ICOL 

IF  (A (ICOL, ICOL) .EQ.O.)  PAUSE  ' Singular  matrix . ' 

PIVINV=1 . /A (ICOL, ICOL) 

A (ICOL, ICOL) =1 . 

DO  16  L=1 , N 

A (ICOL, L) =A (ICOL, L) *PIVINV 

16  CONTINUE 

DO  17  L=1,M 

B (ICOL, L) =B (ICOL, L) *PIVINV 

17  CONTINUE 

DO  21  LL=1,  N 

IF (LL.NE. ICOL) THEN 
DUM=A(LL, ICOL) 

A (LL, ICOL) ~0 . 

DO  18  L=l, N 

A (LL, L) =» (LI , L) -A (ICOL, L) *DUM 

18  CONTINUE 

DO  19  L-l ,  M 

B  (LL, L)  =B  (LL,  L)  -B  (ICOL,  L)  *DUM 

19  CONTINUE 
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ENDIF 

21  CONTINUE 

22  CONTINUE 

DO  24  L=N, 1,-1 

IF (INDXR(L) .NE.INDXC(L) ) THEN 
DO  23  K=1,N 

DUM=A (K, INDXR(L) ) 

A(K,  INDXR(L) ) =A (K, INDXC (L) ) 
A(K, INDXC(L) ) =DUM 

23  CONTINUE 
ENDIF 

24  CONTINUE 
RETURN 
END 


c —  *****★****★**********•★**★**★******★★*★**★*★*■*■★★★★*★★★★*★ 

SUBROUTINE  CTX3D (PLYPRP, CTX) 

c —  formulation  of  three-dimensional  PLY  stiffness  matrix 

REAL  NU21 

DIMENSION  PLYPRP (25)  ,  CTX (6,  6) 

NU21  =  PLYPRP (23) *PLYPRP (22) /PLYPRP (21) 

VE  =  1 . / ( (1 . +PLYPRP (24) ) * (1 . -PLYPRP (24) -2 . *PLYPRP (23) *NU21 ) ) 

c 

CTX (1,1)  =  (1 . -NU21*PLYPRP  (23) ) *VE*PLYPRP (22) 

CTX (1,2)  =  (PLYPRP (24) +NU21*PLYPRP (23) ) *VE*PLYPRP (22) 

CTX (1,3)  =  NU21* (1 . +PLYPRP (24) ) *VE*PLYPRP (21) 

CTX (2,1)  =  CTX (1,2) 

CTX (2,2)  =  CTX (1,1) 

CTX (2, 3)  =  CTX (1,3) 

CTX (3,1)  =  CTX (1,3) 

CTX (3, 2)  =  CTX (2, 3) 

CTX (3, 3)  =  (1 .-PLYPRP (24) **2) *VE*PLYPRP (21) 

CTX (4, 4)  =  PLYPRP (25) 

CTX (5,5)  =  CTX (4,4) 

CTX (6, 6)  =  (1 . -PLYPRP (24 ) -2 . *NU21*PLYPRP (23) ) *VE*PLYPRP (22) /2 . 

RETURN 

END 
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Listing  of  the  material  data  file  "matdat.spr"  for  SPHERE 


1 

T3/5203 [pi/4] 

. 772E11 

. 235E11 

. 27E10 

.  772E11 

. 27E10 

. 13E11 

.  57E10 

. 57E10 

. 268E11 

. 15E10 

. 15E10 

.  40E8 

. 246E9 

.  68E8 

-0.5EO 

- . 80E0 

. 02E-6 

22.5E-6 

.  6E-6 

0.20 

. 181E12 

. 103E11 

.  28E0 

0.52E0 

. 717E10 

2 

H-IM6/EP [pi/4 ] 

. 871E11 

.267E11 

.  44E10 

. 871E11 

.  44E10 

. 154E11 

. 50E10 

. 50E10 

. 302E11 

. 35E10 

. 154E10 

.  56E8 

.  15E9 

.  98E8 

-0.50E0 

- . 80E0 

.0 

0.0 

0.0 

0.1 

. 203E12 

.112E11 

0. 32E0 

0. 60E0 

.  84E10 

3 

E-Glass [pi/4] 

. 223E11 

.  737E10 

.  477E10 

.  223E11 

.  477E10 

. 123E11 

. 319E10 

. 319E10 

. 747E10 

. 106E10 

.  61E9 

.  31E8 

. 118E9 

.  72E8 

-0. 5E0 

80E0 

8.6E-6 

22.1E-6 

0.6E-6 

0.07 

. 386E11 

. 827E10 

0.26E0 

0.60E0 

. 414E10 

4 

A L  2014 

. 101E12 

. 462E11 

. 462E11 

. 101E12 

. 462E11 

. 101E12 

. 276E11 

. 276E11 

. 276E11 

. 449E9 

. 897E9 

.  449E9 

. 897E9 

. 283E9 

-0.5E0 

-0.5E0 

23.2E-6 

23.2E-6 

0.0 

0.0 

. 725E11 

. 725E11 

0.313E0 

0.313E0 

. 276E11 

5 

SS  AM350 

. 284E12 

.  133E12 

. 133E12 

.  284E12 

. 133E12 

. 284E12 

. 75E11 

. 75E11 

. 75E11 

.  114E10 

. 228E10 

. 114E10 

.  228E10 

. 738E9 

-0.5E0 

-0.5E0 

11.7E-6 

11.7E-6 

0.0 

0.0 

. 200E12 

.  20.0E12 

0.318E0 

0.318E0 

.  750E11 

Number  Name 

Cll 

Cl  2 

Cl  3 

C22 

C23 

C33 

C44 

C55 

C66 

X 

X' 

Y 

Y' 

S 

Fxy* 

Fyz* 

alfl 

alf2 

bta2 

Em* 

Ex 

Ey 

nu21 

nu32 

Es 

Cij's  are  the  elements  of  the  effective  stiffness  matrix  in  Pa 

Ei's  are  the  elements  of  the  Young's  Moduli  of  an  unidirectional  ply  in  Pa 
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This  is  the  print  out  of  the  input-output  sequence  for  the  code  "sphere". 
This  run  is  for  a  hybrid  sphere  of  b/a=1.25,  and  50%  of  thickness  of 
T300/5208  inside  and  rest  of  IM6/Epoxy  outside.  The  sphere  is  subjected 
to  internal  pressure. 


INPUT  OUTPUT  FILE  NAME 
T3IM6 

INPUT  INNER  DIAMETER 
7 

>1.0 

INPUT  THE  VALUE  OF  b/a  OF  THE  SPHERE 
>1.25 

INPUT  INTERNAL  PRESSURE  &  EXTERNAL  PRESSURE 

7 

>1  0 

★★♦♦MATERIAL  INFORMATION**** 
INDEX  NUMBER  MATERIAL  DESCRIPTION 


1  T3/5203 [pi/4] 

2  H-IM6/E? [pi/ 4 ] 

3  E-Glass [pi/4] 

4  AL  2014 

5  SS  AM350 

***  INPUT  PLY  LAYUP  INFORMATION  (max  25  layers)*** 

INPUT  NO  OF  LAYERS  (max  25  layers) 

>2 


**  LAYERS  ARE  NUMBERED  FROM  INSIDE  TO  OUTSIDE  ** 

***  FOR  LAYER  1  *** 

INPUT  MATL  INDEX  NO.  &  THICKNESS  FACTOR 
>1  .5 

**  LAYERS  ARE  NUMBERED  FROM  INSIDE  TO  OUTSIDE  ** 

***  FOR  LAYER  2  *** 

INPUT  MATL  INDEX  NO.  &  THICKNESS  FACTOR 
>2  .5 

INPUT  NO  OF  EQUIDISTANT  CALCULATION  POINTS  THRU  THICKNESS  IN  EACH  LAYER 
(max  11,  min  2) 

>3 


***  PLEASE  WAIT  *** 
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LISTING  OF  THE  OUTPUT  FILE  "T3IM6" 

inner  radius (meter) ,  outer  radius (meter) 

.500000  .625000 
b/a  SPHERE 
1.25000 

mat  #  r, meter  sig(r),MPa  sig(th)=sig(phi) ,MPa 
1  0 . 500000E+00  -. 100000E-05  0.244346E-05 

1  0. 531250E+00  -.637914E-06  0.194125E-05 

1  0 . 562500E+00  -.377960E-06  0.162349E-05 

2  0 . 562500E+00  -.377959E-06  0.179410E-05 

2  0. 593750E+00  -.166265E-06  0.159966E-05 

2  0 . 625000E+00  O.OOOOOOE+OO  0.149123E-05 

mat  #  r, meter  eps(r)  eps (th) =eps (phi)  u(r),mm 
1  0 . 500000E+00  -.879822E-10  0.266238E-10  0.133119E-07 

1  0 . 531250E+00  -.577207E-10  0.208252E-10  0.110634E-07 

1  0 . 562500E+00  -.361736E-10  0.170920E-10  0.961423E-08 

2  0 . 562500E+00  -.343097E-10  0.170920E-10  0.961423E-08 

2  0 . 593750E+00  -.192543E-10  0.148012E-10  0.878822E-08 

2  0 . 625000E+00  -.765717E-11  0.134000E-10  0.837503E-08 

mat  #  r (meter)  strength  ratio 
1  0 . 500000E+00  0.240422E+09 

1  0 . 531250E+00  0.307092E+09 

1  0 . 562500E+00  0.348892E+09 

2  0 . 562500E+00  0.327183E+09 

2  0 . 593750E+00  0.368268E+09 

2  0 . 625000E+00  0.272434E+09 

FAILURE  LOCATION  &  FAILURE  PRESSURE 

mat  #  r (meter),  INT  PRESSURE,  EXT  PRESSURE  (in  Pa) 

1  0 . 500000E+00  0.240422E+09  O.OOOOOOE+OO 

END  OF  OUTPUT 


*U,S.  GOVERNMENT  PRINTING  OFFICE:  l*W  -  MR-076/44240 


